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Stability Criteria of 3D Inviscid Shears 

Y. Charles Li 



Abstract. The classical plane Couette flow, plane Poiseuille flow, and pipe 
Poiseuille flow share some universal 3D steady coherent structure in the form of 
£H ■ "streak-roll-critical layer" 1221 |21l 1201 , As the Reynolds number approaches 

infinity, the steady coherent structure approaches a 3D limiting shear of the 
, form (U(y, z), 0, 0) in velocity variables. All such 3D shears are steady states 

I i of the 3D Euler equations. This raises the importance of investigating the 

3 ' stability of such inviscid 3D shears in contrast to the classical Rayleigh theory 

r 1 j ' of inviscid 2D shears. Several general criteria of stability for such inviscid 3D 

^ , shears are derived. In the Appendix, an argument is given to show that a 2D 

■ limiting shear can only be the classical laminar shear. 

• i-H 

>v 

43 . 

i ^~f ' 1. Introduction 

There has been a lot of continuing interest in searching for 3D steady solu- 
J> ■ tions (or traveling wave solutions in a different frame) in plane Couette flow, plane 

Poiseuille flow, and pipe Poiseuille flow [l7][l6][2l][n][22][2Q][6][4][7]. There 
OS . seems to be confirmation of their existence in experiments [7]. Recent numerical 

studies of [22] |21j [20] reveal that the so-called lower branch steady states in the 
plane Couette flow, plane Poiseuille flow, and pipe Poiseuille flow share some uni- 
versal steady coherent structure in the form of "streak-roll-critical layer" . As the 
Reynolds number approaches infinity, the steady coherent structure approaches a 
| 3D limiting shear of the form (U(y, z), 0, 0) in velocity variables. All the 3D shears 

of this form are steady states of the 3D Euler equations. This raises two important 
questions: (1). What is the special property of the limiting shear? (2). What is 
the nature of stability of 3D inviscid shears in contrast to the classical Rayleigh 
theory of 2D inviscid shears? The first question was addressed in [14j . It turns out 
that the limiting shear satisfies a necessary condition: J AU f(U)dydz = for any 
function /. We shall address the second question in this study. We shall use the 
channel flow (plane Couette flow and plane Poiseuille flow) as the example. 

As the Reynolds number decreases from infinity, the limiting 3D shear as a 
steady state deforms into the lower branch steady state; while the 3D shear it- 
self undergoes a slow drifting toward the classical laminar shear. In fact, all the 
shears (3D and 2D) form a stable submanifold of the classical laminar shear. These 
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shears can play a fundamental role in the transition to turbulence from the classical 
laminar shear 



2. Inviscid Channel Flow 

The inviscid channel flow is governed by the 3D Eulcr equations 

(2.1) d t Ui + UjUij = —p t i, u hi = 0; 

where (^1,^2,^3) are the three components of the fluid velocity along (x,y, z) 
directions, and p is the pressure. The boundary condition is the so-called slip 
condition 

(2.2) U2(x, a, z) = 0, U2{x, b, z) = 0; 

where a < b, and ui (i = 1,2,3) are periodic in x and z directions with periods i\ 
and £3. 

We start with the steady shear solutions of the 3D Eulcr equations: 

til = U(y, z), U2 — 0, U3 = 0, p — po (a constant), 

where U (y, z) is periodic in z with period £3. Linearize the 3D Euler equations with 
the notations 

u 1 = U(y,z) + e ik( - x ~ ct) u(y, z) + c.c. , u 2 = e ik( - x ~ ct) v(y, z) + c.c, 



u 3 = e tk{x ~ ct) w(y,z) + c.c, p -> p + e ili{x " ct) p(y, z) + c.c 

where k is a real constant and c is a complex constant, we obtain the linearized 3D 
Euler equations 

(2.3) ik(U ~ c)u + vlly + wll z = -ikp, 

(2.4) ik(U - c)v = -p y , 

(2.5) ik(U - c)w = -p z , 

(2.6) iku + v y + w z — 0. 
Two forms of simplified systems can be derived: 

(2.7) k 2 {U - C)V = Oy [{U - C)(Vy + W Z ) - {UyV + U z w)] , 

(2.8) k 2 {U - c)w = d z \{U - c)(v y +w z )~ (U y v + U z w)} , 

with boundary condition v(a, z) = v(b, z) = and v, w are periodic in z; and 

(2.9) [U-c) 2 V- [(U-c)- 2 Vp] = k 2 p, 

with boundary condition d y p(a, z) — d y p(b, z) — and p is periodic in z. We are not 
successful in utilizing the system (|2.7p - ()2.8j) . System (|2.9p turns out to be fruitful. 
The first result that can be derived from system (|2.9p is the Howard semi-circle 
theorem. 

Theorem 2.1. [8] [3] The unstable eigenvalues (if exist) lie inside the semi-circle 
in the complex plane: 

M + m\ 2 o fM-m\ 2 

+ Cj < , c = c r + iCi, Ci > 0; 



where M = max VlZ U , and m = min a! 2 U. 
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PROOF. Multiply ()2.9|) with p, integrate by parts, and split into real and imag- 
inary parts; we obtain that 

(2.10) / / UGdydz = c r / / Gdydz, 

JO Ja Jo Ja 

(2.11) / / U 2 Gdydz = (c 2 + cf) / / Gdydz, 
where 

Let 
then 



o J a 



G = \U - c\- 4 [\Vp\ 2 + k 2 \p\ 2 } 



M = max U, m — min U, 



/ / (U -m)(M -U)Gdydz>0. 

JO Ja 



10 J a 

Expand this inequality and utilize (|2.10p - (|2.1ip . we arrive at the semi-circle in- 
equality in the theorem. □ 

Our next goal is to find a counterpart of the Rayleigh criterion [2]. For that 
goal, we need to introduce the transform 

(2.12) p={U-c)- x V , 

then p satisfies 

AU 2VJ7-VE/1 



(2-13) ^ + [w^ c -~-W^f\^ k ^ 
with the boundary condition 

(2.14) d y p + jf^P = 0, at y = a, b. 

It turns out that we can only derive results when U satisfies the constraint 

(2.15) U y = 0, at y = a,b; 

in this case, p satisfies the simplified boundary condition 

(2.16) d y p = Q, at y = a, b. 

Theorem 2.2. For U(y, z) satisfying the constraint U y — 0, at y — a, b; if U (y, z) 

has an (inviscid) unstable eigenvalue, then 

(i) 

^(^™)=». 

at some point (y*, z*) in the interior of the domain, for some c (ci > 0) 
in the semi-circle (c r — M±m ) 2 + c 2 < ( M ~ m ) 2 ; where M = max y . z U, 
and m = min y z U; 

(2) 

Id 2 - U 2 

2 17? ^|VC/| 2 + C/AC/>0, 

\U — c\ A 

at some point (y*, z*) in the interior of the domain, for some c (ci > 0) 
in the semi-circle (c r - M±™) 2 + c | < (*Lpi) 2 ; where M = max y , z U, 
and m = min,^ U . 
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Proof. Multiply (|2.13|) with p, integrate by parts, and split into real and 
imaginary parts; we obtain that 



(2.17) 



(2.18) 



U -%AU-2 {U - Cr) 



\U 



u 



-AU 



\U-c 
U-c 



u 



|V[/| ; 



{ p\ dydz — 0. 



l p\ 2 dydz > 0, 



Equation (|2.18[) directly implies the first claim in the theorem. The second claim is 
along the spirit of the Fjortoft theorem [2]. Multiply (j2~T5)) by c r and add (j2~TT)) . 
we obtain the second claim. □ 

Next we will derive a relation between Ci and k. 

Theorem 2.3. For U(y,z) satisfying the constraint U y — 0, at y = a,b; let A = 
max yiZ |V?7|, B = max ajZ |AZ7|; then the unstable eigenvalue (if exists) and the 
wave number k satisfy the condition 

(kcA 2 < 2A 2 + Bc h (Cj>0), 



equivalently, 



kc '-Tk 



Proof. Multiply (|2.13[) with p and integrate by parts, we obtain that 



(2.19) 



J a 



AU 

U-c 



Nu\ 



(u- 



\p\ 2 dydz 



o J a 



[|Vp| 2 + k 2 \p\ 2 ] dydz. 



The left hand side of (|2.19|) is less than or equal to 



B 

W\ 



A 2 



\p\ dydz. 

\Ci\" J Jo J a 

The right hand side of (|2.19p is greater than or equal to 

" ' 1 \p\ 2 dydz. 



Thus 



J a 



B A 2 

h 2^ > fc 2 , (when a > 0); 

Ci cf 



which leads to the claim of the theorem. 



□ 



It is obvious that Theorems 12.21 and 12.31 apply to the 2D shears U(y) too. 
Theorem 12.21 is not the exact 3D counterpart of the 2D Rayleigh criterion. The 
exact counterpart seems elusive. Next we will derive a variation formula for the 
unstable eigenvalue. This type of formulas was initially derived by Tollmien [19] 
[15j for 2D shears. They are useful in deriving unstable eigenvalues near neutral 
eigenvalues. For 3D shears in atmosphere problems pQ, specific approximations 
can make the stability problem very similar to the 2D Rayleigh problem. In such a 
case, a similar variation formula can also be derived to predict unstable eigenvalues 
near neutral eigenvalues pQ. In our current case, no approximation can be made, 
and we have a much harder problem. We have to work with the pressure variable 
of which the singularity nature is not clear even for 2D shears. We can derive a 
formula near an unstable eigenvalue, but its limit to a neutral eigenvalue is elusive 
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and finding a neutral eigenvalue here is more challenging than finding an unstable 
eigenvalue in contrast to the 2D shear problem. Let (p, c, k) and (pi, ci, fci) be two 
unstable eigenfunctions (if exist) to (|2.9|) , then 



(2.20) V- [(U-c)- 2 Vp\ = k 2 (U-c)- 2 p, 

(2.21) V • [(U - ci)- 2 Vpi] = kj(U - ci)- 2 p x , 

Multiply (|2.20[) by pi and (|2.2ip by p, integrate and subtract, we get 



(ci - c) 



o 



ci + c- 2t7 

„ (C/-Ci)2([/-C)2 
6 



[Vp ■ Vpi + k 2 pp\\ dydz 



(fci-fc)(fci + fc) 



/o y, (c^-ci) ; 

from which we obtain the variation formula 



rdydz 



dc 

dk = ~ k 



P 



(2.22) 



J a (U-C) 
s rb 



-dydz 



1 



a (U-CY 



[Vp ■ Vp + k 2 p 2 ] dydz 



The merit of this formula is that it does not involve dp. This formula is valid at an 
unstable eigenvalue (c — c r + ia, Cj > 0). If the unstable eigenvalue lies on a curve 
c = c(k) that leads to a neutral eigenvalue c° = c(k°), = as in the case of a 2D 
shear, then by the semi-circle theorem 12. 11 c° = U (yo, Zq) for some (yo, zo). In such 
a case, the limit k — > fc° of (|2.22[) is still very attractive even though finding the 
neutral eigenvalue c° here is more challenging than finding an unstable eigenvalue in 
contrast to the case of 2D shears. On the other hand, the limit seems very singular 
(even for 2D shears). First of all, in the limit k — > fc°, equation (|2.9|) is singular, so 
the limiting eigenfunction p will be singular too. The following simple equation 

t 2 q" + atq' + (3q — 0, (a, (3 constants) 

already show a variety of singular solutions near t — 0. Moreover, the singularities 
generated by (U — c°) in the integrals in (|2.22|) add to the challenge. 

Similarly, starting from (|2.13|) . (|2.15p and (|2.16|) , we can derive the following 



AU 



(ci - c) 
= (ki-k)(h + k) 

from which we obtain the variation formula 



2VC/ • VC/(ci + c - 2U) 



(U-c)(U-c) 

b 

ppxdydz 



(U -a) 2 (U -c) 2 



ppidydz 



dc 
dk 



2k 



(2.23) 



7 

■J a 


) 

p 2 dydz 




f 


' AU 


4VU ■ VZ7 


J a 


_{U-c) 2 


(U - c) 3 



p 2 dydz 
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3. Viscous Channel Flow 

The viscous channel flow is governed by the Navier-Stokes equations 

(3.1) d t Ui + UjUi^ = —p ti + euijj, Ui t i = 0; 

where (^1,^2,^3) are the three components of the fluid velocity along (x,y,z) 
directions, p is the pressure, and e = 1/R is the inverse of the Reynolds number R. 
The boundary condition is 

(3.2) u\ (x, a, z) = a, u±(x, 6, z) — (3, Uj(x,a,z)=Uj(x,b,z)=0,(J = 2,3); 

where a < b, a < (3, and Ui (i — 1,2,3) are periodic in x and z directions with 
periods i\ and £3. For the viscous channel flow, the 3D shears mentioned above are 
no longer fixed points, instead they drift slowly in time (sometimes called quasi- 
steady solutions): 

(e £ ' A [/(y,z),0,0). 

By ignoring the slow drift and pretending they are still fixed points (or by using 
artificial body forces to stop the drifting), their unstable eigenvalues will lead to 
transient nonlinear growths as shown numerically [12j. The corresponding linear 
Navier-Stokes operator at (U(y, z), 0,0) is given by the following counterpart of 



(3.3) ik(U - c)u + vU y + wU z = -ikp + e[A - k 2 ]u, 

(3.4) ik(U - c)v = -p y + e[A - k 2 ]v, 

(3.5) ik(U - c)w = -p z + e[A - k 2 }w, 

(3.6) iku + v y + w z =0. 

Again two forms of simplified systems can be derived: 

d y {[A-k 2 - ikR(U - c)] (v y +w z ) + ikR{vU y + wU z )} 

(3.7) ~k 2 (A-k 2 )v + ik 3 R{U -c)v = 0, 

d z {[A ~ k 2 - ikR(U - c)] (v y +w z )+ ikR(vU y + wU z )} 

(3.8) -k 2 (A-k 2 )w + ik 3 R(U - c)w = 0, 

with boundary condition 

v(a, z) = v(b, z) = v y (a, z) = v y (b, z) — w(a, z) — w(b, z) = 0, 
and v,w are periodic in z] and the other form 

(U - c) 2 V • {{U - c)" 1 [e(A - k 2 ) - ik(U - c)] -1 Vp} = ikp 

(3.9) +4(A - fc 2 )V • {[e(A - k 2 ) - ik{U - c)] _1 Vp} , 

the boundary condition of which is complicated. 

The system (|3.7|l - (|3.8|l looks quite convenient for numerical simulations. Mul- 
tiply (|3.7p by v and ()3.8|) by w, integrate and add the two equations, we obtain the 
following expression for the eigenvalue in term of the eigenfunction: 

-[A + E + H + kRRe(G)}, c r = —[B + Im(G)], 



kRD L v ' D' 
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where 



A = J [\Vv y \ 2 + \\7w z \ 2 + k 2 \v y \ 2 + k 2 \w z \ 2 ] > 0, 

B = ju{\ Vy \ 2 + \ Wz \ 2 + ,y Wz+Vy ^ 

+k 2 \v\ 2 + k 2 \w\ 2 ) , real; 



D = J \ \Vy\ 2 + \W Z \ 2 + VyW Z + VyW Z 

+k 2 \v\ 2 + k 2 \w\ 2 ^j > 0, 
E = J [vyyW Z y + v yz w zz + k 2 v y w z + c.c] , real; 
G = —i J (v y + w z )(vU y + wU z ), complex; 
H = k 2 I \\Vv\ 2 + |Vw| 2 l +fc 4 / \\v\ 2 + \w\ 2 } > 0. 



Notice that the expression of c r has no explicit dependence upon the Reynolds 
numbber R, but it does depends on R implicitly via the eigenfunction. An unstable 
eigenvalue corresponds to kci > 0. Without loss of generality, assume k > 0. Notice 
that A + E > 0, and 

\G\<9 J[\v y \ 2 + \w z \ 2 + \v\ 2 + \w\ 2 }, 

where 

(3.10) g = ma^{\\U y \\ L o ,\\U z \\ L o a }. 

Then from the expression of Ci , one can obtain the following theorem by comparing 
the terms H and kR Re(G). 

Theorem 3.1. Let g be given by &3.10\) , (k, c) be an eigenmode; if Rg < min(fc, k 3 ), 
(k > 0); then Ci < 0, i.e. a stable eigenvalue. 

The type of claims in Theorem 13.11 and their improvements have been investi- 
gated intensively for 2D shears 18 [9\ |10j . On the other hand, as mentioned at the 
beginning, we are more interested in unstable eigenvalues and for large Reynolds 
number R as in 



4. Appendix: 2D Limiting Shear Can Only Be the Classical Laminar 

Shear 

We take the plane Couette flow as the example, for plane/pipe Poiseuille flow, 
the argument is the same. Assume that as the Reynolds number R — > +oo, a 
2D steady state of the plane Couette flow approaches the limiting shear (U(y), 0). 
Denote by 

(U(y) + u(x,y),v(x,y)) 
the steady state which is periodic in x. Assume that (it, v) and their spatial deriva- 
tives are of order o{R~ 1 ^ 2 ) as R — > +oo. Then to the leading order 0(R^ 1 ), 
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ud x u + vd y u = -d xP + -dp, 

Ud x v = -d v p, 
d x u + d y v = 0. 

Taking an average in x— direction (over the period), we get 

vdyU = ^d 2 y U, 
= -dyf, 

d y v = 0. 

By the boundary condition of v in y direction, TJ = 0. Thus 

d 2 y U = 0, i.e. U = ci + c 2 y. 

That is, U has to be the laminar linear shear. 

As shown in |14j . the corresponding 3D limiting shear U(y, z) does not have to 
be the linear shear, rather satisfies a constraint 

J AUf(U) dydz = 

for any /. In fact, the lower branch 3D limiting shear is far away from the linear 
shear [2 2) [21] |20j . Does this hint that the lower branch steady state does not 
exist in 2D? Numerical simulations could not find any 2D steady state other than 
the linear shear [5]. 
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